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ABSTRACT 

SNe la are good distance indicators because the shape of their light curves, which can 
be measured independently of distance, varies smoothly with luminosity. This suggests that 
SNe la are a single family of events. Similar correlations are observed between luminosity 
and spectral properties. In particular, the ratio of the strengths of the Si II A5972 and A6355 
lines, known as 7vL(Si II), was suggested as a potential luminosity indicator. Here, the physical 
reasons for the observed correlation are investigated. A Monte-Carlo code is used to con- 
struct a sequence of synthetic spectra resembling those of SNe with different luminosities 
near B maximum. The influence of abundances and of ionisation and excitation conditions 
on the synthetic spectral features is investigated. The ratio 7\L(Si II) depends essentially on the 
strength of Si II A5972, because Si II A6355 is saturated. In less luminous objects, Si II A5972 
is stronger because of a rapidly increasing Si Il/Si III ratio. Thus, the correlation between 
1Z{$>\ II) and luminosity is the effect of ionisation balance. The Si II A5972 line itself may 
be the best spectroscopic luminosity indicator for SNe la, but all indicators discussed show 
scatter which may be related to abundance distributions. 
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1 INTRODUCTION 

Type la supernovae (SNe la) form a relatively homogeneous su- 
pernova subclass. They supposedly originate from accreting C-0 
white dwarfs that undergo explosive burning as they appro ach the 
Chandrasekhar mass limit dHillebrandt & Niemeverj|200d) . Their 
enormous light output is powered by radioactivity of synthesised 
56 Ni and of its daughter nuclide 56 Co dArnettlll982h . 

Differences in the synthesised 56 Ni mass result in a spread i n 
luminos i ty of ~3 mag in the B band ( iMazzali et aT1l200lL [2007h . 
IPhillipsN 19931) found that the decline in B-band magnitude within 
15 days after B maximum, Amis (B), correlates inversely with the 
B-band maximum luminosity M s. max and thus can be used as a lu- 
minosity indicator. Thus, SNe la can be used as distance indicators 
in cosmological studies ("standard candles"). 

Optical spectra reflect in detail the variations among SNe la. 
In the photospheric phase, which lasts until a few weeks after 
maximum light, P-Cygni features (blueshifted absorption plus less 
prominent emission centred at th e rest wavelength) are superim- 
posed on a "pseudo-continuum" jPauldrach et al.|[l996t) . The line 
positions and strengths depend on the abundances a nd on the ex- 
citation and ionisation conditions within the debris. iNugent et al.l 



dl995h showed that the ratio of the depths of two Si II lines (A5972 
and A6355) at B-band maximum, 72.(Si II), is larger for dimmer 
SNe, and suggested that this is due to temperature differences. Ac- 
cordingly, 72.(Si II) could be used as a distance-independe nt spec- 
troscopic luminosity indicator (e.g. iBongard et al.ll200ol) . Other 
spectral features ha ve subsequently been suggested as po ssible lu- 
minosity indicators dHachinger. Mazzali & BenettiirJoOol he reafter 
IPaper j IBongard" et al.ll2006klFolev, Filippenko & Jhall2008h . Still, 
the detailed physics behind the 72.(Si II)-L correlation have not been 
explained. 

This is the purpose of this work, where we explain the 
7Z(Si II)-L correlation investigating the details of the connection 
between the physics of the envelope and line strengths. We present 
a sequence of spectral models calculated with a Monte-Carlo (MC) 
radiative transfer code. The synthetic spectra resemble real SNe 
with different luminosities shortly after B max. In the model at- 
mospheres, we analyse the conditions which influence the line 
strengths entering 72.(Si II). 

In the following, first we describe the concept of the synthetic 
spectral sequence (Sec. 2). In Sec. 3, we show that our models re- 
produce the TZ(Si ll)-L correlation. In Sec. 4 we explain the correla- 
tion analysing the details of the formation of the lines contributing. 
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After discussing spectral luminosity indicators and other applica- 
tions of line-strength measurements (Sec. 5), we give conclusions 
(Sec. 6). 



2 CONCEPT OF THE MODEL SEQUENCE 

2.1 Model SN envelope and radiative transfer simulation 

Our wor k makes use of a ID Monte Carlo (MC) rad i ative trans- 
fer code dAbbott & Lucvl[l98l iMazzali & Lucvlll993L iLucvll 19991 
and lMazzalill2000l . which describes the version used in this work) 
to compute synthetic SN spectra. The code has succe ssfully been 
applied to many SNe (SN 1990N: lMazzali et al.lll993l etc.). 

The code computes the radiative transfer through the SN 
above an assumed p hotosphere. The envelope is assumed to expand 
homologously (e.g. iRopke & Hillebrandlll2005l) . meaning that the 
ejecta move radially with r = v • t, where r is the distance from 
the centre, t the time from explosion and v the velocity. Radius and 
velocity can therefore be used interchangeably for each gas parti- 
cle. The density structure is taken from a s tandard explosion model 
(W7. lNomoto. Thielemann & Yo koi 1984), scaled according to the 
time t. The chemical composition of the model atmosphere can be 
freely adjusted, but is assumed to be homogeneous in radius. 

From the photosphere, which is located at an adjustable u p h, 
"packets" of thermal radiation [JJ" = B V {T P \,)] are emitted into 
the atmosphere. As they propagate through the atmosphere, radia- 
tion packets can undergo Thomson scattering on electrons or line 
excitation-deexcitation processes, which are treated in the Sobolev 
approximation. The process of photon branching is included, which 
implies that the transitions for excitation and deexcitation can be 
different. In all processes, radiative equilibrium is conserved. Tak- 
ing the backscattering into account, the code matches a given output 
luminosity L. 

The action of the radiation field onto the gas is calculated us- 
ing a modified nebular approximation, which mimics the effects of 
NLTE. The envelope is discretised into 80 shells. For each shell, a 
radiation temperature Tr and a dilution factor W are calculated. 
These quantities mostly determine the excitation and ionisation 
state. 

The code iterates the radiation field and the gas conditions. 
After convergence, the output spectrum is obtained from a formal 
integral solution of the transfer equation dLucvlll999l) . 

2.2 Sequence of synthetic spectra 

For our study we computed a sequence of 25 spectral models. 
The start and end points of the sequence were chosen to match 
a luminous SN la (model no. 1: SN 2003du, Ami 6 (B)= 1.02, 
MB.modd = —19.16) and a dim one (model no. 25: SN 1999by, 
Ami5(B)= 1.90, MB >m0 dci = —16.86), respectively, close to B 
maximum. Observed and synthetic spectra are shown in Fig.[H fur- 
ther data and references are given in Table 1 . 

The free model parameters shown in Tableland Fig.|2]were 
iterated to optimise the fit. Some deviations between models and 
data remain, but they have no influence on our investigation. In 
the luminous models there is an excess flux in the red caused by 
the assumption of a blackbody lower boundary. Dim models show 
a mismatch in line velocities, which sh ould vanish with a refine d 
density and abundance stratification (cf. lTaubenberger et alfcooih . 

The focus of our fits is on obtaining a reasonable temperature 
structure. The decrease of temperature from luminous to dim SNe 
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Figure 1. Selected spectra from the model sequenc e, and observed counter- 
parts (black) with Ami5(£?) values (see lPaperB) . All spectra are plotted 
logarithmically and have been adapted to the total reddening of SN 2003du. 
Approximate identifications are given for features relevant for this work. 



Table 1. First (no. 1) and last (no. 25) model of the sequence: Data of the 
objects and spectra used for modelling. 

Object/epoch/ Amis Total reddening Distance modulus Redshift References 
model no. [mag] E(B — V)^ [mag] [mag] z forSN,/x,2 



03du/+1.2d/#l 1.02 0.01 32.79 0.0068 Stanishev et al. (2007) 

99by/+3.9d/#25 1.90 0.01 30.75 0.0021 Garnavich et al. (2004) b) 

a * The values given represent the galactic extinction E(B — V)g^, as the host galaxy extinction 

is approximately zero for both SNe ( see references). E{B — V)u U \: from NED (see 

acknowledgements); we averaged the Burstein & Hciles i 1982) and Schlegel, Finkbeiner & Davis 
1 19981) values. 

b * Spectrum combined from original spectra with different grating setups (taken between JD 
2451312.64 and 2451312.66, available at http://www.cfa.harvard.edu/supernova/) and 
fine-calibrated to match photometry by multiplying with a linear function (using IRAF and 
SYNPHOT routines/data, see acknowledgements). 



not only implies that the spectra appear redder, but also influences 
line strengths and 72.(Si II) via excitation and ionisation. Choos- 
ing appropriate models for SNe 03du and 99by as the ends of the 
sequence, we ensure that we explore a realistic range of tempera- 
tures. The temperature structure of our models can to some degree 
be influenced by the choice of photospheric velocity. We fitted the 
pseudo-continuum of both SNe as well as possible. We did not in- 
crease temperatures further for SN 03du, because this would make 
Si III A4560 too strong. The resulting effective temperatures (cf. 
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Table 2. Overview of the final parameters for models no. 1 and 25. Mass fractions less significant have been omitted. 

Model lg I ^ J i?ph t a Element abundances (mass fractions) 

^ @ ' [kmr 1 ] [d] X(OJ X(Mg) X(Si) X(S) X(Ca) X(Ti) X(Cr) X(Fe) b) X( 56 Ni) b) 

no. l:03du/+1.2 9.61 8900 21.0 0.056 0.13 0.56 0.19 0.0075 0.0006 0.0006 0.0075 0.060 

no. 25: 99by/+3.8 8.82 6500 21.0 0.68 0.085 0.21 0.016 0.0008 0.0016 0.0016 0.0002 0.0025 

a ' Time from explosion onset. The rise times for 03dti and 99by are assumed to be 19.8d and 17. Id, respectively. This complies with observational studies i Riess et al. 1999: Con lev et al. 2006) and leads to a 
similar time from explosion of ~21 .0d for both models. 

b ' The abundances of Fe, Co and Ni in our models are assumed to be the sum of 5 ^Ni and its decay chain products (^ 6 Co and 5 ^Fe) on the one hand, and directly synthesised / progenitor Fe on the other hand. 
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Figure 2. Abundances of the 99by and 03 du model. The Fe and Ni abun- 
dances given are X(Fe)o and X( 56 Ni)o (see tab. [2}- T ne SN 99by model is 
generally characterised by very small abundances of burning products, ex- 
cept for relatively large amounts of Ti and Cr (Fig. [5). Despite the decrease 
in Si abundance, Si II A5972 is stronger in dimmer objects. 



order to match the evolution of the S II "W-trough" better, it de- 
creases with a larger slope among luminous SNe and with a smaller 
slope among dim ones. 

We assign Amis {B) values to the models no. 1 and 25 and to 
four intermediate ones. In order crudely to assign Amis (B) values 
to intermediate models, we created "optimal-fit" models of the SNe 
2002bo [Ami 5 (B)= 1.17], 1994D [Ami 5 (B)= 1.32], 1991M 
[Ami 5 (B)= 1.51] and 1993H [Ami 5 (B)= 1.70] as reference 
models. We then compared the sequence models to these by calcu- 
lating a weighted sum of quadratic deviations in the input param- 
eters (Lboi, v P h and abundances), where the deviation in Lboi was 
given the strongest weight. The Ami5(B) values of the reference 
models were then assigned to the four sequence models that gave 
the closest match. 

Fig. [D shows the spectral models assigned a Ami^jB) 
value with their o bserved counterpa rts (SN 2002bo: iBenetti et al .1 
|2004 SN 1994D: IPatat et alj[l996l SN 1991M: Asiago archive, 
SN 1993H: CTIO archive, with kind permission of M. M. Phillips). 
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Figure 3. Input parameters L^ \ and t> pn of our models, and resulting ef- 
fective temperature T e ff. The photospheric blackbody temperature Tph after 
convergence is given for comparison. 
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(<tb: Stefan-Boltzmann's constant; t ~ 21. Od: time from explo- 
sion) of our models are shown in Fig. [3] For comparison we also 
give the photospheric temperatures T p h which are calculated by the 
code, taking the backscattering into account. 

The other models of the sequence (no. 2-24) were calculated 
decreasing log L in constant steps. This leads to a relatively uni- 
form distribution of models in Mb ■ The photospheric velocity and 
all element abundances were changed by constant amounts from 
one model to the next. An exception is the sulphur abundance: in 



3 R(SI II) IN THE MODEL SEQUENCE 
3.1 Line-strength measurements 

In order to compare models and observations, we measured the 
strengths of Si II A5972 and Si II A6355, and 7£(Si II) in synthetic 
and observed spectra. 

To quantify the line strength, the flux level inside a feature is 
compared to an estimated "pseudo-continuum", which we define 
as a linear function approximately tangential to the spectrum at the 
edges of the feature. We then measure the line strength as fractional 
depth: 



FD(Xi 



Fpc(Xfd) - F{X FD ) 
Fpc(Xfd) 



(1) 



where F(Xfd) and Fpc(Xfd) are the actual flux and the 
pseudo-continuum, respectively, at a wavelength Xfd, which is 
chosen so as to obtain the maximum To avoid the influence 
of noise, measurements were made in spectra convolved with a nor- 
malised Gaussian function (FWHM: 15A). 



3.2 Evolution of FDs and ft(Si II) 

Fig. [4] shows our measurements of FD(Si II A5972) and 
FD(Si II A6355). While FD(Si II A5972) clearly increases 



1 iNugentetafU 19951) quantified line strength as "depth", i.e. the absolute 
difference between the pseudo-continuum and the actual flux. The original 
7Z(Si II) depth ratio of Si II AA5972, 6355 behaves like the corresponding 
FD ratio in observations. 
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Figure 4. A5972 and A6355 FD in models and observed reference spectra. 
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Figure 5. Ratio of Si 11 A5972 to Si 11 A6355 FD [~7£(Si II)]. 

with decreasing luminosity, FD(Si II A6355) is approx- 
imately constant in all measurements. The resulting ratio 
FD(Si II A5972)/FD(Si II A6355) ~ K(Si n) is shown in Fig. [5] 
The value increases towards lower luminosities. 

Observed spectra are less well behaved than the synthetic 
ones, and show some constant deviation at the dim end. This could 
be improved wit h more sophistic ated models using abundance 
stratification (e.g. IStehle et alJ2005n . possibly including some con- 
tribution of Na I D to the A5972 feature in dim SNe [See e.g. 
iFilippenko et al .1 dl992l) . Na I would make too broad a contribution 
to the A5972 feature in our homogeneous models.]. Yet, the qual- 
itative evolution of the two Si II features, as well as their different 
behaviour, is correctly reproduced. 



4 ANALYSIS 

Here we analyse the behaviour of 7Z(Si II), which is rooted in the 
evolution of Si II A5972 and Si II A6355. We show that the latter 



line is easily saturated because it is very strong. Si II A5972, in 
contrast, is stronger in dimmer objects owing to a larger fraction of 
Si II. 

Furthermore, we in vestigate if Ti II lines h ave an influence on 
72.(Si II), as proposed bv lGarnavich et all d2004l) . 

4.1 Quantifying the conditions for absorption in our models 

Absorption strengths are determined by intrinsic line strengths, 
abundances, density, ionisation and excitation structure. For ab- 
sorptions which are (mostly) caused by a single ion, line formation 
can be conveniently analysed. Before performing the analysis, we 
introduce our method and quantities important in this context. 

Region of the envelope relevant for absorption at Xfd 

The reference wavelength of a FD, Xfd, is blueshifted with respect 
to the rest frame wavelength A rcs i, because absorption occurs in ma- 
terial moving towards the observer. Therefore, as seen from the ob- 
server's point of view, the processes relevant for a particular FD 
take place within a thin layer where the atoms have a line-of-sight 
(ray-parallel) velocity i)| = (1 — ^ £C -) ■ c. In the narrow-line limit 
(Sobolev approximation), with v <x r, this zone is an infinitely thin 
"iso-D| | plane" perpendicular to the ray. 

The code calculates physical quantities describing absorption 
in each radial shell. Thus, a weighted average over several shells 
is required to obtain values corresponding to a measured FD (and 
the respective iso-Un plane). This is denoted by angle brackets "< 
... >" below. The weight of each shell is determined by the number 
of photons impinging on the absorption surfaces of the multiplet 
within the shell. 



Sobolev optical depth and absorption fraction 

Both of the features in 7£(Si II) are doublets (AA5958, 5979 and 
AA6347, 6371, respectively). The total fraction of photons ab- 
sorbed ("absorption fraction" 21), which is expected to correspond 
to FD values, can then be calculated as: 

»:= <i-II e ~ T! ' fc >' (2) 

I 

where ri_k is the Sobolev optical depth for line I in shell k, 
defined as 

^=*-^i^(i-£L)- (3) 

Here, n, and rij are the number density of ions in the lower 
and upper level of the line, i and j, respectively, and t is the time 
from explosion. All other symbols have their usual meaning. 

The influence of ionisation and occupation numbers 

Besides the intrinsic line strength Bij, we identify two major in- 
fluences on Ti,k\ First, the number density of Si II, iVsiii, which 
depends on the ionisation balance. Second, the level occupation 
numbers, which determine how many of the ions are available for 
the transition. The influence of occupation numbers can be encap- 
sulated into an "excitation factor" (symbols as above): 
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Figure 6. Absorption conditions for Si II A5972. 

Upper panel: Ratio of Si II /Si III ionisation fractions averaged over the 
shells. 

Middle panel: (p(Si II)}, excitation factor Si and the resulting op- 
tical depth (t);. (t); is roughly determined by (p(Si II)) and Si 
(crudely (r); oc (p(Si II)) ■ Si), (t); and Si are plotted for one represen- 
tative line (A5958); the data for the other doublet line are approximately 
proportional. 

Lower panel: Absorption fraction 21, and measured FD. Measurements in 
this and all following plots have been performed several times. Fluctuations 
in FD and also in 21 (due to fluctuations in A p u ) are represented by error 
bars (whenever large enough). 



"Sin \h V 9i n i ]k J 

This contains the fraction of ions in the line's lower level and 
the correction for stimulated emission, which usually varies rela- 
tively weakly. 
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Figure 7. Absorption conditions for Si II A6355; plot analogous to Fig. [6] 
The range of absorption fractions and FDs is small because of saturation. 
Apart from intrinsic line strengths, the physics behind Si II A5972 and Si II 
A6355 is similar. Minor differences depend on different blueshifts of the 
features (see Appendix[g] Fig.|B5j: Si II A5972 forms at >9000 km s" 1 
(luminous models) to ~8000 km s (dim models); the stronger A6355 
line always forms well above the photosphere at ~12500 km s _1 . 

~90% towards lower luminosities and temperatures. Yet, the opti- 
cal depth of Si II A5972 increases (filled squares / dotted line). 

The absorption fraction 21 (lower panel, open squares / dark, 
dotted line) increases towards lower luminosities, following the op- 
tical depth, as 21 = (1 — Yli e~ n - h ) ~ Ti.t) for small optical 
depths (first order). Measured line FDs (circles / grey line) reflect 
this trend. 

The offset between the exact values of FD and the absorp- 
tion fraction, which increases somewhat towards the dim end, is 
discussed in Appendix lAl As it turns out, ~50% of this offset is 
due to overlap with the emission part of the P-Cygni feature. The 
strength of this generally correlates with the absorption. 



4.2 Si II A5972 - the trend behind the H(S\ \i)-L correlation 

The Si II A5972 absorption is stronger in dimmer SNe. This be- 
haviour, which is the basis for the 7£(Si II) sequence, is due to the 
Si II to Si III ratio (Fig. [6] upper panel), which increases by a factor 
of ~ 100 from luminous to dim SNe. Other ionisation stages are 
negligible in the layers probed. 

The middle panel shows in detail the influence of the increas- 
ing Si II density (crosses / dark, solid line) and of the excitation con- 
ditions on the optical depth. The fraction of ions in the lower level 
of the line (~ excitation factor, crosses / dashed line) decreases by 

© 2008 The Authors. Journal compilation © 2008 RAS, MNRAS 



4.3 Si II A6355 - a line with constant strength 

Si II A5972 evolves with luminosity, while Si II A6355 does not 
(Fig. [7] bottom panel). This is due to saturation: because Si II 
A6355 has large optical depths (middle panel), the absorption frac- 
tion (bottom panel) approaches unity, and it varies only weakly 
(21 = (1 — Yli e~ T '- k ) ~ 1 for large r). Measured FDs show a 
shallow trend, as do absorption fractions. The FD values indicate, 
however, an additional amount of residual flux. Again, we show in 
Appendix |A] that 50% of this amount is due to the P-Cygni emis- 
sion. 
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Model no. 25 
Modified model -X(Ti)*10 




Figure 8. Model no. 25 and a modified version with 10-fold Ti abundance, 
where Ti lines contribute in the order of 10% to the strength of the 5972A 
feature. A multitude of Ti II transitions is active, with strongest lines at 
AA5888.89, 5910.05, 5987.39, 5994.94, 6001.4, 6012.75. FD(Si II A5972) 
does not yet show an increase, also because Si II A5972 gets weaker due to 
the slightly modified temperature structure. The larger Ti abundance leads, 
however, to an unrealistic opacity increase below ~5200A. 



4.4 Alternative explanations for FD (Si II A5972) and 

ft(Si II)? 

iGarnavich et all d2004 suggested that the evolution of Si II A5972 
is due to a contribution of Ti II lines in dim SNe. However, this 
possibility can be ruled out reviewing the list of active lines. In 
order for Ti II lines to make a significant contribution, the Ti abun- 
dances must increase by a factor of > 10 in our SN 1999by model. 
This however leads to exceedingly large strengths for many other 
Ti lines (Fig. [8}. We also do not find any further species that have 
a significant influence on the evolution of Si II AA5972, 6355, with 
the possible exception of Na I in dim SNe (see Sec. 13. 2b . 



5 DISCUSSION 

5.1 Spectroscopic luminosity indicators: 1Z(Si II) and Si II 
A5972 line strengths 

Having clarified the behaviour of observed FD values, we now 
briefly discuss the use of line strengths and ratios as luminosity 
indicators. We have shown that FD(Si II A5972) correlates with 
luminosity and drives TZ(Si II). Thus, FD(Si II A 5972) it self (or 
similarly the equivalent width EW, as measured in lPaper 3) should 
be the most direct and meaningful luminosity indicator. 

There are, however, uncertainties even in FD(Si II A5972) 
as a spectroscopic luminosity indicator. First, the trend in the line 
is due to the combination of trends that are nonlinear in tempera- 
ture, luminosity and also in Amis(B) in the relevant ranges (see 
Sec. |4). In order to obtain luminosities from spectra, one has to 
choose an appropriate functional form (e.g. linear or second order) 
to fit the relation between feature strength and luminosity. Inaccu- 
racies introduced by the use of such a functional form can easily 
be overlooked or misinterpreted as scatter. Second, intrinsic vari- 
ations in the SNe and their spectra limit the precision of spectro- 
scopic luminosity indicators. Differences in the abundance distri - 
butions among equally luminous SNe [see e.g. lTanaka et al.l l f2008l) . 
iMazzali & Podsiadlowskil d2006l) l are a possible reason for this (see 
also below, Sec.l5.2t. 



One way to reduce the error would be to use a number of spec- 
troscopic luminosity indicators. Even quantities which are theoreti- 
cally more elusive than FD may be useful, because different indica- 
tors respond differently to spectral peculiarities. Besides ratios like 
TZ(Si II), EW is also worth considering, because it accounts better 
for absorption in the line wings. SNe for which spectroscopic lu- 
minosity determination is unreliable could be identified if different 
indicators gave inconsistent luminosities. 

The accuracy of spectroscopic luminosity determination could 
be improved if additional luminosity indicators were found that 
rely on species other th an Si II. A number of such quantities 
has already been defined l Nugent et al ] | 19951 , iBongard et al]|2006| , 
Paocr I iFolev. Filippenko & Jhal 120081) . However, the validity of 
these potential luminosity indicators still needs to be critically as- 
sessed in terms o f the physical differ ences among SNe with differ- 
ent luminosities jMazzali et alj20o"7h . 



5.2 FD(Si II A5972) and SN physics 

The Si II A5972 optical depths, absorption fractions and FDs are 
quite tightly correlated. The line thus reflects the physics in the Si- 
dominated layers of most SNe la much better than Si II A6355. It 
would be interesting to compare the Si II A5972 profiles in B max. 
spectra of several SNe with similar luminosity and explore the dif- 
ferences. If the Si distribution peaks in the zone with the largest 
Si II ionisation fraction, one may expect relatively large FDs. On 
the other hand, a smeared-out Si distribution, or one with several 
peaks mig ht cause an unusually shallow feature, a s observed in 
SN 2006X llWang et al.l200lj,lElias-Rosa et alj|2008l (in prep.)]. 

Going one step further, it would be interesting to infer the 
distribution of Si using st ratified abundance models ( Steh le et al .1 
l2005t I Mazzali et al. 1 120081) . Th is would be especia lly valuable in 
the context of SN la theory (see lMazzali et al.ll2007h and explosion 
physics. 



5.3 Other features evolving with luminosity 

Apart from Si II A5972 and A6355, other features in photospheric 
spectra of SNe la seem to evolve with l uminosity. The Si II A 4130 
line behaves similar to Si II A5972 l lBronder et al .1 1200 8h . and 
the S II trough is strongest for SNe with intermediate luminosity 
jPaper j) . In Appendix|B] these two features are investigated. 



6 CONCLUSIONS 

Based on a sequence of 25 synthetic spectra constructed to fit SNe 
la with different luminosities, we have been able to explain the 
7£(Si n) - L relation. 

The increasing strength of Si II A5972 is responsible for the in- 
crease of 7?,(Si II) from luminous to dim objects. This is in turn due 
to ionisation: in dimmer SNe, Si II A5972 is stronger because the 
Si II fraction increases more strongly than than the level occupation 
numbers decline. In contrast, Si II A6355 is saturated because of its 
large optical depths. Thus, the line has a rather constant strength 
and only a minor influence on 7£(Si II). 

An investigation as performed here shows how observed line 
strengths can be related to SN physics. Si II A5972 probably pro- 
vides the most reliable information about the Si layer in optical 
spectra. 

If luminosity is to be inferred from spectra, a simultaneous 
evaluation of several line strengths or ratios is most promising. In 
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order to find additional reliable luminosity indicators, one can ap- 
ply our methods to further lines caused by single ions. 
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Figure Al. Contributions to the residual flux and continuum definitions. 
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APPENDIX A: THE OFFSET BETWEEN FRACTIONAL 
DEPTHS AND ABSORPTION FRACTIONS 

While the general trends in Si II A5972 and A6355, and thus in 
1Z(S\ II), can be explained by absorption effects, the mismatch be- 
tween 21 and FD deserves closer investigation. Here, we analyse 
this offset. 

Al Influences on the fractional depth besides absorption 

There are two potential reasons for a difference between absorption 
fraction and FD: First, photons may be shifted into the observed 
wavelength of a feature in re-emission processes, and enhance the 
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residual flux. Second, the estimated pseudo-continuum may devi- 
ate from the flux behind the absorption surfaces, as seen from the 
observer. 

To trace these effects, we counted photon packets represent- 
ing different components of the residual flux in our simulations 
and recorded the pseudo-continuum estimate for each FD measure- 
ment. 



Components of the residual flux 

The residual flux at a measurement wavelength \fd can be split 
into contributions of different origin (Fig. lAll l, 

One contribution ("penetrating flux" F p ) directly depends on 
the absorption fraction 21. It consists of photons crossing the mul- 
tiplet's absorption surfaces without being absorbed, and is approx- 
imately given by 



F P ^{1~%)F EC 



(Al) 



where Fec ("effective continuum") is the flux impinging on the ab- 
sorption surfaces of the multiplet, multiplied by a correction factor. 
This factor, usually a bit smaller than unity, accounts for the fact 
that photons passing the absorption surface do not have a 100% 
chance of escaping the envelope. 

The rest of the residual flux ("emission" F t ) consists of pho- 
tons which have been shifted to an observer-frame wavelength 
^\fd at the absorption surfaces (corresponding to \fd) or even 
closer to the observer. In the following, photons that underwent ab- 
sorption and reemission at the absorption surfaces will be called 
"P-Cygni emission" {F et pc yg ). The amount of P-Cygni emission 
will generally correlate with the absorption strength. The flux of 
photons that were reprocessed into an observer-frame wavelength 
~Afc (for example in line processes) at depths closer to the ob- 
server will be called "additional reprocessed flux" (F eia d). It is an 
upper limit estimate for the flux that enters the feature without cor- 
relation to the absorption conditions. 



The relation between 21 and FD 

In general, the relation between absorption fraction, emission and 

1 - FD, is: 



measured residual flux -J^- 

Fpc 



F PC 



(1-FD): 



F p + F c 
Fpc 



rpc 



Fpc 



(A2) 



2t and FD are closely related only if jpf*- and the additional 
reprocessed flux are roughly constant with luminosity. As the data 
will show, this is the case for Si II A5972 and A6355, but not always 
for the lines discussed in Appendix iBl 



A2 Si II A5972 

For Si II A5972, emission causes a moderate offset between pen- 
etrating flux and measured residual flux F res = 1 — FD (Fig. lA2l 
squares / dark line and circles / grey line, respectively). On average, 
about half of the emission (triangles / dark, dotted line) is P-Cygni 
emission, which increases with decreasing luminosity. The addi- 
tional reprocessed flux, in turn, is relatively constant (not explicitly 
shown), as is F fi^ L . Thus, the measured FDs correlate well with the 
absorption fractions. 

Still, the offset between FDs and absorption fractions has a 
practical consequence: It makes the relative decrease in FD from 
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Figure A2. Contributions to the emission, and FD for Si II A5972. All 
fluxes are given relative to the pseudo-continuum. The penetrating flux 
(1 — 21) X Z, EC and the emitted flux J'' can be added up on the left 

^ PC PC 

scale to obtain a prediction for the residual flux (not explicitly shown). This 
coincides roughly with the measured residual flux (grey curve). The grey 
curve represents the fractional depth with units on the right hand side. 
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Figure A3. Emission contributions and FD for Si II A6355; plot analogous 
to Fig.lA2l 



intermediate to very luminous SNe larger, increasing the influence 
of Si II A5972 on ft(Si II). 



A3 Si II A6355 

Intuitively, the residual flux in the saturated Si II A6355 line should 
be lower than observed. Again, this is in roughly equal parts due to 
P-Cygni emission and the additional reprocessed flux. Both these 
contributions and are practically constant within the luminos- 
ity sequence. 



A4 Concluding remarks 



For the two lines determining 7?.(Si II), jfi^- is constantly « 1 and 
the additional reprocessed flux is small (<10% of the residual flux). 
Even under slightly different simulation conditions, we expect FDs 
to correlate well with the absorption in their respective lines, be- 
cause redwards of Si II A5972, the number of lines which might 
cause deviations is relatively small. 
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Figure Bl. Absorption conditions for Si II A4130 (for log(L bo i/LQ) > 
9.1, see text). 

Upper panel: Ionisation and abundance trends: Si mass fraction, and Si II 
and Si III ionisation fractions (shell-averaged). 

Lower panel: (p(Si II)} and excitation factor £j; (t); and absorption frac- 
tion 21. (t); and £; are plotted for one representative line (A4128); the data 
for the other triplet lines are approximately proportional. 

APPENDIX B: INVESTIGATION OF FURTHER 
FEATURES (SI II 4130A, S II 5640A) 

SN la spectra contain a wealth of features, of some of w hich show 
evolution with luminosity (see e.g. lNugent et al.l 19951 and lPaperll) . 
Calcium absorptions frequently seem to suffer from saturation 
and high-velocity contrib utions jHatanoetalJll999l :l Mazzali et al.l 
l2005l : iTanaka et alj|2006h . Interesting trends have been found for 
Si II A4130 and the S II "W-tr ough": Si II A4130 correlates 
with luminosity like Si II A5972 dBronder et~ai]|2008l) . although 
the Si II A4130-lum inosity correlation suffers from more scatter 
dBronder et alj|2008l . Fig. 10). The S II trough, in contrast, shows 
a parabo lic trend, with largest strength at intermediate luminosities 
dPaper J) . Both features are largely caused by a single ion and can 
be investigated with the technique applied to 7£(Si II). 



Bl Si 11 A4130: a supplementary luminosity indicator 

The data for this feature are summarised in Figs. lBll and lB2l The 
line can be measured only in spectra with log(L/ Lq) >9.1. At 
lower luminosities, the Ti II lines at ~4000-430fjA typical of dim 
SNe affect the feature. 

Among normal SNe, the absorption fraction shows a trend 
similar to Si II A5972 (Fig. |BT1 lower panel). Measured FD val- 
ues (Fig. IB2I circles) correlate with the absorption fractions, but 
increase only slowly towards lower luminosities. This is mainly 
due to an increase in the additional reprocessed flux (not explic- 
itly shown), which leads to an increase in the emission (triangles 
/ dark, dotted line). Contamination effects (e.g. by S II A4162.7, 
Ni II A4067.0 and Co II A4160.7 ) are a possible reason for this. 

As noted by iBronder et all d2008h . the Si II A4130 strength 
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Figure B2. Contributions to the emission, and measured residual flux (grey 
curve) for Si II A4130. The grey curve represents the fractional depth with 
units on the right hand side. Measurements of the residual flux in the ob- 
served reference spectra are shown for comparison. 

could be used as a luminosity indicator for intermediate to bright 
objects. However, FD does not follow the absorption fraction very 
well for the feature. Thus, the relation b etween luminosity and Si II 
A4130 line strength shows large scatter llBronder et al.l d2008h . Fig. 
10], and we regard the feature as a supplementary luminosity indi- 
cator at most. 

B2 S II A5640: a mix of temperature and abundance trends 

The blue part of the S II W-trough has differen t depths in 
the sy nthetic and observed spectra (Fig. \T\ see also iKotak et"aT] 
l2005h . Thus, we chose to investigate only the red component 
(lines/multiplets AA5606, 5640). The pseudo-continuum was set as 
if the whole trough were measured. 

The results are shown in Figs. lB3l and lB4l Sulphur has differ- 
ent level structures and a lower ionisation threshold than silicon, 
so that the S II/S III ratio is always > 0.5. Furthermore, because 
of the steeply declining line velocity (Fig. lB5b . the temperatures in 
the layers probed do not vary as much as for the Si lines within the 
sequence. 

In the range from large to intermediate SN luminosity, a mild 
increase in the S II fraction (Fig. IB3I upper panel, asterisks / dark, 
dashed line) increases the optical depths somewhat. At low lumi- 
nosities, the S II fraction approaches 100%, so that decreasing S 
abundances (Fig. IB3I upper panel, crosshairs / solid line) and level 
occupation numbers lead to a decreasing absorption strength (Fig. 
IB3I lower panel). 

The resulting FDs (Fig. IB4I circles) in the synthetic spectra 
increase from large to intermediate lu minositie s and then decline 
again, similar to the parabolic trend in lPaper 3 for the EW of the 
W-trough. The decline in FD from intermediate- to low-luminosity 
SNe is amplified by an increase in the additional reprocessed flux, 
which affects dim SNe more strongly (triangles / dark, dotted line). 
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Figure B5. Line-of-sight velocities corresponding to \fd values (in the 
synthetic spectra) for all features discussed in this paper. Weaker features 
tend to form near the photosphere at lower velocities, and trace the photo- 
spheric recession (filled squares / black, dashed line). 



Figure B3. Absorption conditions for S II A5640 (red part of W-trough). 
Abrupt variations in the values are due to jumps in the measurement posi- 
tions (see Fig. lB5t . 

Upper panel: S mass fraction of our models, and S II and S III ionisation 
fractions (shell-averaged). 

Lower panel: (p(S II)) and excitation factor (r); and absorption frac- 
tion 21. (t)i and Si are plotted for one example line (lower level 13.7eV, 
non-metastable). The lines of the S II A5640 multiplet originate from dif- 
ferent levels at ~13.7...14eV, including a metastable one. Trends in (t); 
values for other lines will thus deviate somewhat, causing the absorption 
fraction to deviate from the individual optical depth trends. 
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Figure B4. Emission contributions and FD for S II A5640. Plot analogous 
to Fig.lBll 
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